#ifndef __POTENTIALS__
#define __POTENTIALS__

#include "stdio.h"
#include "stdlib.h"
#include "complex.h"
#include "math.h"

#define PI (3.1415926535897932)

/**
 * The following four functions are
 * vector and scalar potentials 
 * that comes from the literatures
 */

void f_xcL0(double complex *fxcL0, double *rho_0, int Nx, int Ny, int Nz);

void f_xcLT(double complex *f_xcL, double complex *f_xcT, 
	    double *rho_0, double complex omega_input, 
	    int iVK, int Nx, int Ny, int Nz);

void f_xc_ALDA(double complex *f_xc, 
	       double *rho_0, double *QD_epsilon, 
	       int Nx, int Ny, int Nz);

void eta_xc_xi_xc(double complex *eta_xc, double complex *xi_xc, 
		  double complex *fxcL0, double complex *f_xcL, double complex *f_xcT, 
		  double *rho_0, double complex omega_input, int Nx, int Ny, int Nz);

#endif
